Von Karman Vortex Street within an Impacting Drop 
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The splashing of a drop impacting onto a liquid pool produces a range of different sized micro- 
droplets. At high impact velocities, the most significant source of these droplets is a thin liquid jet 
emerging at the start of the impact from the neck that connects the drop to the pool. We use ultra- 
high-speed video imaging in combination with high-resolution numerical simulations to show how 
the ejecta gives way to irregular splashing. At higher Reynolds number, its base becomes unstable, 
shedding vortex rings into the liquid from the free surface in an axisymmetric von Karman vortex 
street, thus breaking the ejecta sheet as it forms. 

PACS numbers: 47.55. D-, 47.32. ck, 47.20. Ma, 47.20. Dr 
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Liquid drop splashing is part of our daily lives, from 
the morning shower to natural rain Q. While it 
has been studied for more than one hundred years [3|, 
it is only recently that advances in high-speed imaging 
techniques 0, Q have revealed its early dynamics Q Q . 
Splashing refers herein to the breakup of a drop into 
smaller droplets during impact. Understanding the un- 
derlying mechanism which produces the smallest droplets 
is important for example for the number of microscopic 
aerosols which remain when those satellite droplets evap- 
orate. Such aerosols affect human health and can act as 
nucleation sites during cloud formation. 

For high speed drop impact on a liquid pool, the ejecta 
sheet is the first stage leading to splashing. It was first 
observed in the inviscid numerical simulations of Weiss 
& Yarin [8] and in the experiments of Thoroddsen 
When the drop impacts at higher velocity, the speed of 
these ejecta sheets increases and they become thinner. 
The radial stretching of the sheets reduces their thickness 
even further and they can remain intact even at thick- 
nesses well under a micron Q. When they eventually 
rupture, they can produce a myriad of very fine spray 
droplets. However, this mechanism does not continue for 
ever, at a critical Reynolds number the smooth ejecta 
gives way to a more random splashing, which counterin- 
tuitively may produce fewer small droplets. 

To understand the mechanisms leading from contin- 
uous ejecta sheets to irregular splashing, a systematic 
study of the early dynamics was conducted with ultra- 
high-speed video imaging, over a range of impact veloci- 
ties [/, liquid viscosities /i and droplet diameters D [lot . 
Figure [T] shows a classification of the results in terms of 
Reynolds number Re = pDU / ji, where p is the liquid 
density, and splashing parameter which relates to the 
Weber number We = pDU'^/a, where a is the surface 
tension, i.e. K = Wey/Re. We are interested here in the 




FIG. 1. Characterization of the ejecta regimes. (#) Smooth 
ejecta sheet. (■) Irregular splashing, (v) Bumping, (a) 
Quartering. (♦) Protrusions rising up along the side of the 
drop. 



higher K regime, where splashing occurs 

The classification in Fig. [T] focuses on the ejecta 
shapes. In the lower range of Re (more viscous liquids), 
a smooth ejecta sheet emerges between the drop and the 
pool (#). However in the highest range of Re, isolated 
droplets emerge from the neck, followed by a disturbed 
hquid surface, and no coherent ejecta can be identified, 
i.e. irregular splashing occurs (■). In the intermediate 
regime (Re ^ 2000 — 6000), the ejecta sheets show a large 
variety of repeat able shapes. We have grouped them into 
3 classes. At lower K (lower impact velocities, ), surface 
tension prevents the formation of an ejecta sheet. How- 
ever, we observe some protrusions travelling up along the 
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FIG. 2. Comparison between experiment and axisymmetric 
numerical simulation for a bumping case. U = 4.04m/s, Re = 
3.55 X 10^ ,K = 7.44 x 10"^. From top to bottom, observations 
at time 30, 80, 130, 180 and 230 /is after contact, (a) Experi- 
mental observation. The static dark points correspond to dust 
on the camera sensor. The video was taken at 200,000 frames 
per second, (b) Numerical simulation of the drop impact for 
the same conditions as in the experiment presented in (a). In 
the last image, the leading part of the eject a sheet becomes 
smaller than the grid size, by stretching between the apex and 
the tip, and thus breaks into non-physical droplets. The ax- 
isymmetric simulations cannot include the three-dimensional 
effects, such as the breakup of the tip observed in (a). Sup- 
plemental videos show the two evolutions ^IQ] . 



side of the drop, without ejection of droplets outwards. 
At higher K (a)^ the ejecta sheet is more developed; how- 
ever, it stays attached to the drop, stretching the ejecting 
sheet between the expanding tip of ejecta and the drop 
entering the pool. This stretching can lead to the explo- 
sive rupturing of the sheet, which generates fast droplets 
of a large range of sizes through slingshot [6] . In the up- 
per range of (v) we observe an intriguing phenomenon 
where the free-standing sheet interacts strongly with the 
downwards moving drop surface. This is shown in the se- 
quence of Fig. [2ja), herein referred to as the bumping of 
the ejecta. The ejecta is strongly bent by the drop, and 
then folds at its apex. Overall snapshot of a bumping 
ejecta was included in Thoroddsen et at. [5] [their Fig. 
8(c)]. 

Those experimental results clearly show the effect of 
the Reynolds number on the transition towards irregular 
splashing. Moreover, it shows new dynamics of the ejecta 
sheet interacting with the drop. It suggests that those 
interactions could underlie the irregular splashing. 

To test this idea, we have chosen to reproduce the im- 
pact by numerical simulations. It is only recently that nu- 



merical simulations managed to identify the ejecta sheet 
d, Q, [Hi, due to the extreme range of scales involved, 
and the challenges of interfacial flow simulations (l6j . 
The intricate shapes observed herein were thus not seen 
numerically until now. We use the freely available code 
Gerris 17H19|. for its high parallelization and dynamic 
adaptive grid refinement, which allow us for the first time 
to reach enough precision to fully resolve the dynamics 
of the ejecta. This code uses the volume- of -fluid method 
to solve the incompressible Navier- Stokes equations. 

Axisymmetric simulations faithfully reproduced all of 
the experimentally observed features, as we demonstrate 
in Fig. [2] for the bumping case. The shape of the drop in 
the simulation is perfectly spherical, ruling out the hy- 
pothesis that small deviations from spherical drop shapes 
in the experiments could be responsible for the drop in- 
teraction with the ejecta sheet. 

To study this transition to irregular splashing we in- 
crease i?e, while keeping K constant, from a smooth 
ejecta sheet [Fig. [H#)] to irregular splashing (■). This 
was done for two different K values, corresponding to the 
bumping (v) and quartering (a) regimes. 

The position of the base of the ejecta rx [Fig. [Ha,b)] 
follows very closely the geometric relation predicted by 
Josserand & Zaleski [TBI, rx = C'rj, independent of i?e, 
where r j is the radius where an unperturbed drop would 
meet the pool surface. Pure kinematics suggests that 
the ejecta emerges at an angle intermediate between the 
drop surface and the pool, thus should increase with rj. 
A simple geometric model Q suggested that increases 
as ^ ~ where t* is the time nondimensionalized by 
r = D/U ^ whereas the simulations show that grows 
linearly before bumping [Fig. [3fc)]. However, the ejecta 
rises faster for higher Re. The collapse of the curves 
in Fig. [3jd) shows that grows at a rate proportional 
to \/Re. The angle of the ejection velocity vector at 
the middle of the base also follows a similar trend, and 
increases proportionally to ^/Re. 

At lower i?e, 6 increases slowly enough for the ejecta 
to escape the drop. However, from Re > 3000, the ejecta 
sheet rises too fast, thus impacting the drop surface. The 
resulting bumping sharply decreases 0. This interaction 
of the drop and the ejecta sheet observed experimentally 
occurs earlier at higher Reynolds numbers, eventually 
breaking the ejecta sheet. This is consistent with the 
interpretation that this interaction is responsible for the 
irregular splashing observed at higher Re. 

Vorticity also appears to play an important role in 
the dynamics of the ejecta sheet [Fig. [3fe,f)]. The flow 
around a curved free surface generates vorticity (see for 
instance Batchelor [20] 5.14). Therefore, concentrated 
vorticity is observed in the simulations near points T and 
B at the base of the ejecta, as the flow moves around 
the base to enter the ejecta. At the early stage of the 
ejecta formation, both sides of the base produce a sim- 
ilar strength of vorticity. This initial vorticity scales as 
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FIG. 3. Evolution of the base of the eject a sheet with Re at 
K = 7.44 X 10^. Quantities are nondimentionahzed by the 
drop diameter D, the drop impact velocity U and the drop 
entry time r = D/U . (a) Definition sketch. The base of the 
eject a sheet is defined as the segment between the two points 
of maximum curvature of the interface (T on the drop side and 
B on the side of the pool). The angle of the ejecta sheet is 
the angle between the horizontal and the normal to the base. 

(b) Evolution of the ejecta base radial position, defined as the 
distance from the axis of symmetry to the middle of TB in 
(a), vs. the nondimensional time t* = t/r, for Re from 1,000 
to 6,000. The solid curve is 1.23rj, where rj — '\/t*(l — t*) 
is the radius where an undisturbed drop would meet the pool. 

(c) Evolution of (in degrees) for Re from 1,000 to 6,000. (d) 
Same curves as (c), where the angle is scaled by \/Re. (e,f) 
Evolution of the maximum positive vorticity (red) and maxi- 
mum absolute negative vorticity (blue) in the liquid near the 
ejecta base for Re = l, 000 (e) and Re = 4, 000 (f). The posi- 
tive maximum is located near T, and the negative maximum 
near B. 



\/Re as observed previously [li]. However, the differ- 
ence in vorticity (absolute values) between the two sides 
is observed to increase initially linearly with time, before 
decreasing again. Moreover, this difference is higher for 
larger Re [Fig. [3];f)]. 

By looking closely at the neck region during the im- 
pact, we can identify fundamental changes in the vortic- 
ity structure as Re is increased (Fig. |4]). Note that in Fig. 
m^a) most of the liquid in the sheet originates from the 
pool, in agreement with dye visualizations [9]. For the 
lower range of Re [Fig. |31^a,d)], the vorticity stays con- 
centrated near the free surface at the neck of the ejecta 



sheet. As there is stronger vorticity generated at the top 
of the ejecta base, a vorticity layer of one sign separates 
the drop and the pool liquids but it remains stable. K 
affects the shape of the outer part of the ejecta sheet, as 
we observe by comparing Fig. |4l^a,d), consistently with 
experimental observations [6|. For intermediate Re^ the 
interface remains stable in its early evolution. In the 
bumping case (b), the rising ejecta sheet contacts the 
downwards moving drop surface. This creates a shear 
instability, generating a toroidal vortex structure around 
the entrapped bubble. In the quartering case (e), the 
ejecta sheet leaves the neck region to climb up the drop, 
pulled by higher surface tension. This also creates a shear 
instability between the climbing liquid from the pool and 
the drop liquid moving down, forming a row of vortex 
rings of the same sign. These vortices near the free sur- 
face leave their signature [2l| by creating waves below 
the rising sheet, a feature also observed experimentally. 
However, all such vortical effects are absent from inviscid 
theory and simulations [1, [22| 

At even higher Re [Fig. |4l^c,f)], vorticity is shed behind 
the base of the ejecta sheet, in a way reminiscent of the 
von Karman vortex street, here forming alternating-sign 
vortex rings. For the first 7 shedding cycles, the local 
Reynolds number based on the radial speed and width of 
the neck takes value around 70 and the Strouhal number 
St = fD/U is around 0.11±0.05, in good agreement with 
related Karman streets. Moreover, the free surface at 
the base is also affected by the shedding, reinforcing the 
oscillations [23]. During the early shedding [Fig. IHg)], 
surface tension effects are higher because of the sharper 
surface geometry. As the angle of the neck increases, 
the amplitude of the oscillations increases. The ejecta 
can then climb on the drop at lower K (f), or impact 
alternatively on the drop and the pool (c), in a similar 
_to the bumping, entrapping a row of bubble rings 
Four bubble rings can be clearly identified in Fig. 
m^c), and a fifth one being created. Only well resolved 
bubbles and droplets (larger than 30 cells) are kept in 
the numerics, suggesting that smaller bubbles could be 
entrapped earlier. 

From systematic experimental observations, repro- 
duced with axisymmetric simulations, we have detailed 
a new mechanism explaining the irregular splashing of 
a water drop. Previously studied mechanisms have de- 
scribed the droplet separation from the rim of the ejecta 
BSisI, or the destabilization of a liquid sheet (HUsl. 
Our mechanism however explains the breakup of the 
ejecta sheet by the destabilization of its base, through 
vortex shedding from the free surface. 
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FIG. 4. Vorticity structures during drop impact near the transition regime. In the top images we differentiate the hquids 
originating from the drop (red) and from the pool (blue) , from the air (green) , as can be done in experiments by seeding one or 
the other with fluorescent dye p. The bottom images show the corresponding vorticity in the liquid, (a-c) Bumping transition, 
with K = 7.44 x 10^, for increasing Reynolds numbers: Re = 1000, 3552 & 14500 respectively, (d-f) Quartering transition, 
with i^^ = 3 X 10"^, for Re = 1000, 3552 k 10000 respectively. To allow direct comparison, the images of the first row correspond 
to the same nondimensional time as the ones in the second row, with the same field of view, (c) and (f) correspond to a water 
drop impacting at 2.84 m/s and 1.98 m/s respectively, (g) Details of the early vortex shedding in the same case as (c). The 
period of this shedding shown here is approximately 2 /xs, over a radial distance of 50 /xm. 
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SUPPLEMENTAL MATERIAL 



Experimental setup: The liquid viscosity was var- 
ied in our experiments by using glycerin/ water mixtures 
of various mass fractions, 0, 40, 50, 60, 65, 70, 75 and 
80% of glycerin. The impact velocity was changed by 
releasing the drops from between 0.16-2.16 m height. By 
varying the circular steel nozzle size, we used four differ- 
ent drop diameters. Further details of the experiments 
can be found in Thoroddsen et al. [6]. 

Numerical simulations: Direct comparison between 
the experimental observations and numerical results with 
Gerris [l9| were done on a set of different impact condi- 
tions. A convergence study of a bumping case showed 
that a minimum cell size 1,000 times smaller than the 
drop diameter was necessary to observe the bumping 
shapes, explaining why previous studies did not observe 
this evolution. However, a much finer grid is needed to 



reproduce the evolution of very thin ejecta sheets. This 
was accomplished using up to 16 levels of refinement, 
which means separating the domain in two smaller cells 
16 times in each direction (the size of the smallest cell 
is 2^^ times smaller than the domain). In our simula- 
tions, this corresponds to 28,800 cells per drop diameter. 
At this level of refinement, an equivalent uniform grid 
would have more than 4 billion cells ((2^^) ), while our 
simulation have of the order of 20 million cells. 

At the highest level of refinement, the ejecta sheet can 
become as thin as 500 nm (only 3 cells) near the tip, due 
to the extreme stretching. This value is consistent with 
indirect experimental measurements 

We present here results at level 14, on 64 processors, al- 
lowing a systematic investigation of the parameter space. 

The curvature of the tracer isolines was estimated by 
fitting a Bezier curve of order 3 through 5 points. For 
smaller curvature cases, the number of points was in- 
creased to 11, improving the estimate of the position of 
maximum curvature. 



